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ABSTRACT 

Many proposals have already been made for realizing pro¬ 
grammable matter, ranging from shape-changing molecules, 
DNA tiles, and synthetic cells to reconfigurable modular 
robotics. Envisioning systems of nano-sensors devices, we 
are particularly interested in programmable matter consist¬ 
ing of systems of simple computational elements, called par¬ 
ticles , that can establish and release bonds and can ac¬ 
tively move in a self-organized way, and in shape forma¬ 
tion problems relevant for programmable matter in those 
self-organizing particle systems (SOPS). In this paper, we 
present a general algorithmic framework for shape forma¬ 
tion problems in SOPS, and show direct applications of this 
framework to the problems of having the particle system 
self-organize to form a hexagonal or triangular shape. Our 
algorithms utilize only local control, require only constant- 
size memory particles, and are asymptotically optimal both 
in terms of the total number of movements needed to reach 
the desired shape configuration. 

1. INTRODUCTION 

Imagine that we had a piece of matter that can change 
its physical properties like shape, density, conductivity, or 
color in a programmable fashion based on either user input 
or autonomous sensing. This is the vision behind what is 
commonly known as programmable matter. Programmable 
matter has been the subject of many recent novel distributed 
computing proposals, ranging from shape-changing molecules, 
DNA tiles, and synthetic cells to reconfigurable modular 
robotics. Each of these proposals pursued solutions for spe¬ 
cific application scenarios with their own, special capabilities 
and constraints. 
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We envision systems of nano-sensors devices that will have 
very limited computational capabilities individually, but which 
can collaborate to reach a lot more as a collective. Ideally, 
those nano-sensor devices will be able to self-organize in or¬ 
der to achieve a desired collective goal without the need 
of central control or external (in particular, human) inter¬ 
vention. For example, one could envision using a system 
of self-organizing nano-sensor devices to identify and coat 
(and possibly repair) leaks on a nuclear reactor without 
the need for human intervention; self-organizing systems of 
nano-sensor devices could also be used to monitor environ¬ 
mental and structural conditions in abandoned mines, on the 
exterior of an airplane or spacecraft, bridges and other struc¬ 
tures, possibly also self-repairing the structure— i.e., real¬ 
izing what has been coined as ’’smart paint”. The applica¬ 
tions in the health arena are also endless, e.g., self-organizing 
nano-sensor devices could be used within our bodies to de¬ 
tect and coat an area where internal bleeding occurs, elimi¬ 
nating the need of immediate surgery, or they could be used 
to identify and isolate tumor/malignouos cells. In many ap¬ 
plications, there may be a specific shape that one would like 
the system to assume (e.g., a disc, or a line, or even any 
compact shape). 

Hence, from an algorithmic point-of-view, we are inter¬ 
ested in programmable matter consisting of systems of sim¬ 
ple computational elements, called particles , that can estab¬ 
lish and release (communication or physical) bonds and can 
actively move in a self-organized way, and in general shape 
formation problems in those self-organizing particle systems 
(SOPS). 

1.1 Geometric Amoebot model 

We will use the geometric amoebot model presented in [7| 

[6] as our basic model for SOPS. 

In all of our shape formation algorithms, the set of par¬ 
ticles will maintain a connected structure at all times. We 
assume that we have a graph G(V, E ) that represents the 
relative positions that a connected set of particles may as¬ 
sume — i.e., V represents all possible positions of a par¬ 
ticle (relative to the other particles in their structure) and 
E represents all possible transitions between nodes. In the 
geometric amoebot model we assume that G = G e q t, where 


G e qt is the infinite regular triangular grid graphQsee Part 
(a) of Figure [l]). 

We briefly recall the main properties of the geometric 
amoebot model. Each particle occupies either a single node 
or a pair of adjacent nodes in G eq t, and every node can be 
occupied by at most one particle. Two particles occupying 
adjacent nodes are connected , and we refer to such particles 
as neighbors. 

Particles move through expansions and contractions : If a 
particle occupies one node (i.e., it is contracted ), it can ex¬ 
pand to an unoccupied adjacent node to occupy two nodes. 
If a particle occupies two nodes (i.e., it is expanded ), it can 
contract to one of these nodes to occupy only a single node. 
Performing movements via expansions and contractions may 
represent the way particles physically move, or may be seen 
as a logical ’’look-ahead and then move” logical operation. It 
has several advantages, including allowing particles to abort 
a movement if there is a conflict (see [t] for more details). A 
particle always knows whether it is contracted or expanded 

— in the latter, it also knows along which edge it expands 

— and this information will be available to neighboring par¬ 
ticles. A handover allows particles to stay connected as they 
move; two scenarios are possible: a) a contracted particle p 
can ’’push” a neighboring expanded particle q and expand 
into the neighboring node previously occupied by q , forc¬ 
ing q to contract, or b) an expanded particle p can ’’pull” 
a neighboring contracted particle q to a cell occupied by it 
thereby expanding that particle to that cell, which allows p 
to contract to its other cell. In part(b) of Figure [l] we il¬ 
lustrate a set of particles (some contracted, some expanded) 
on the underlying graph G eq t- 

Particles are anonymous but the bonds of each parti¬ 
cle have unique labels, which implies that a particle can 
uniquely identify each of its outgoing edges. Moreover, for 
each particle the bonds are labeled in a consecutive way in 
clockwise direction so that every particle has the same sense 
of clockwise direction, but the particles may not have a com¬ 
mon sense of orientation in a sense that they have different 
offsets of the labelings (see Figure [l] Part (c)). Each particle 
has a constant-size local memory in which it can store some 
bounded amount of information, and any pair of connected 
particles has a bounded shared memory that can be read 
and written by both of them and that can be accessed using 
the edge label associated with that connection. We assume 
the standard asynchronous model from distributed comput¬ 
ing, where the system of particles progresses by perform¬ 
ing atomic actions, each of which affects the configuration 
of one or two particles. Whenever a particle is activated 
(i.e., performs an atomic action), it can perform an arbi¬ 
trary bounded amount of computation (involving its local 
memory as well as the shared memories with its neighbor¬ 
ing particles) followed by no or a single movement. A round 
is over once every particle has been activated at least once. 

1.2 Our Contributions 

In this paper, we present a general algorithmic framework 
for shape formation problems in SOPS, which constitutes of 
two basic algorithmic primitives: the spanning forest primi¬ 
tive and the snake formation primitive. We present concrete 
applications of these two primitives to two specific shape 
formation problems, namely to the problems of having the 


system of particles self-organize to form a hexagonal shape 
and to form a triangular shape. Both the hexagonal shape 
and the triangular shape formation algorithms are optimal 
with respect to work , which we measure by the total num¬ 
ber of particle movements needed to reach the desired shape 
configuration, as we prove in Theorems [l] and [2] Our algo¬ 
rithms rely only on local information (e.g., particles do not 
have ids, nor do they know n, the total number of particles, 
or have any sort of global coordinate/orientation system), 
and require only constant-size memory particles. 


1.3 Related Work 

Many approaches related to programmable matter have 
recently been proposed. One can distinguish between active 
and passive systems. In passive systems (e.g., DNA com¬ 
puting 11 [2 16 , tile self-assembly systems d[T3][T7]),) the 
particles either do not have any intelligence at all (but just 
move and bond based on their structural properties or due 
to chemical interactions with the environment), or they have 
limited computational capabilities but cannot control their 
movements. We will not describe passive models in detail 
as they are only of little relevance for our approach. On 
the other hand in active systems , computational particles 
can control the way they act and move in order to solve a 
specific task. Robotic swarms, and modular robotic systems 
are some examples of active programmable matter systems. 

In the area of swarm robotics it is usually assumed that 
there is a collection of autonomous robots that have lim¬ 
ited sensing, and communication ranges, and that can freely 
move in a given area. They follow a variety of goals, includ¬ 
ing for example shape formation problems (e.g., [9j [l4] ). 
Surveys of recent results in swarm robotics can be found 
in | 111[T2]. While the analytical techniques developed in the 
area of swarm robotics and natural swarms are of some rel¬ 
evance for this work, the individual units in those systems 
have more powerful communication and processing capabil¬ 
ities than in the systems we consider. 

The field of modular self-reconfigurable robotic systems 
focuses on intra-robotic aspects such as the design, fabri¬ 
cation, motion planning, and control of autonomous kine¬ 
matic machines with variable morphology (see e.g., 10, 19]). 
Metamorphic robots form a subclass of self-reconfigurable 
robots that share some of the characteristics of our geo¬ 
metric model [5]. The hardware development in the field 
of self-reconfigurable robotics has been complemented by a 
number of algorithmic advances (e.g., j3j[l5, 14|), but so far 
mechanisms that automatically scale from a few to hundreds 
or thousands of individual units are still under investigation, 
and no rigorous theoretical foundation is available yet. 

The nubot model |18[ |4] aims at providing the theoretical 
framework that would allow for a more rigorous algorith¬ 
mic study of biomolecular-inspired systems, more specifi¬ 
cally of self-assembly systems with active molecular com¬ 
ponents. While bio-molecular inspired systems share many 
similarities with our SOPS, there are many differences — 
e.g., there is always an arbitrarily large supply of ’’extra” 
particles that can be added to the system as needed, and 
the system allows for an additional (non-local) notion of 
rigid-body movement. 


2. SHAPE FORMATION 


1 The triangular grid graph G eqt is the dual graph of a regular In this paper we focus on solving shape formation problems 

hexagonal tiling in 2D space. in the geometric amoebot model starting from any initial 
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Figure 1: Part (a) shows a section of G eq t; nodes of G e q t are shown as black circles. Part (b) shows five particles on 
G eq t: the underlying graph G eq t as a gray mesh; a particle occupying a single node is depicted as a black circle, and 
a particle occupying two nodes is depicted as two black circles connected by an edge. Part (c) depicts two particles 
occupying two non-adjacent positions on G e qt; the particles have different offsets for their head bond labelings. Part 
(d) shows an intermediate configuration of the HEX algorithm. The seed is depicted in green, retired particles are 
black, and roots are red. Particle p is the last added particle to the retired structure. Hence, edge i connects p to the 
retired particle p' (edge i has the flag p'.snakedir). The red arrow depicts the process of setting p.snakedir in clockwise 
manner for p. 


connected configuration of particles. We present a general 
algorithmic framework for shape formation problems and 
then specifically we investigate the Hexagonal Shape For¬ 
mation (HEX) and the Triangular Shape Formation (TRI) 
problems where the desired shape is a hexagon and a tri¬ 
angle respectively. We formally define a shape formation 
problem as a tuple M — (X, (J) where X and Q are sets of 
connected configurations. We say X is the set of possible 
initial configurations and Q is the set of goal configurations. 
Accordingly, for the HEX problem, Q would be all config¬ 
urations where the positions of the set of particles induce 
a hexagon on G e qt (note that depending on the number of 
particles the constructed shape may not necessarily be a 
perfect hexagon since the outer layer of the hexagon may 
not be fully complete). Similarly, for the TRI problem, Q 
is equal to the set of all configurations that constitute a 
triangle in G eq t (except for possibly the outer layer of the 
triangle, which may be partially full). We say an algorithm 
A solves a shape formation problem M if for any execution 
of A on a system in an arbitrary configuration from X, A 
terminates (i.e., the execution eventually reaches a config¬ 
uration in which each particle does not move anymore) in 
one of the valid configurations in Q. 

Before we proceed, we provide some preliminaries. For all 
algorithms we assume that there is a specific particle we call 
the seed particle, which provides the starting point for con¬ 
structing the respective shape. If a seed is not available, one 
can be chosen using the leader election algorithm proposed 
in [7] . We define the set of states that a particle can be in 
as inactive , follower , root , and retired. Initially, all particles 
are inactive, except the seed particle, which is always in a 
retired state. In addition to its state, each particle p may 
maintain a constant number of flags in its shared memory. 
For an expanded particle, we denote the node the particle 
last expanded into as the head of the particle and call the 
other occupied node its tail. In our algorithm, we assume 
that every time a particle contracts, it contracts out of its 
tail. Note that with this convention, the node occupied by 
the head of a particle still is occupied by that particle after 
a contraction. Part (c) of Figure [l] shows an example of the 
labeling of the heads of two particles on G eq t- 

Generally speaking, the shape formation algorithms we 
propose for hexagonal and triangular shapes progress as fol¬ 
lows. Particles organize themselves into a spanning set of 


disjoint trees where the roots of the trees are non-retired 
particles adjacent to the partially constructed shape struc¬ 
ture (consisting of all retired particles). Root particles lead 
the way by moving in a predefined direction around the cur¬ 
rent structure. 

The remaining particles (i.e., the followers) follow behind 
the leading root particles, hence the system flattens out to¬ 
wards the direction of movement. Once the leading particles 
reach a valid position where the shape can be extended (fol¬ 
lowing the rules for the snake formation for the particular 
shape), they stop moving and change their state to retired 
as well. This process continues until all particles become 
retired. Note that the spanning forest component of this 
general approach is the same for any shape formation al¬ 
gorithm: It is only in the rules that determine the next 
valid position to be filled in the shape structure being built 
that the respective algorithms differ. We determine the next 
valid position to be filled sequentially following a snake (i.e., 
a line of consecutive positions in G eq t), that fills in the space 
of the respective shape structure and scales naturally with 
the number of particles in the system. 

2.1 Spanning Forest Algorithm 

The Spanning Forest algorithm primitive, given in Algo¬ 
rithm [l] is a building block we use for all of our shape forma¬ 
tion problems. This primitive was also used in 7], where we 
present a preliminary self-organizing algorithm for forming 
a straight line of particles. We present the algorithm here 
for completeness. Each particle p continuously runs Algo¬ 
rithm [l] until it becomes retired. If particle p is a follower, 
it stores a flag p.parent in its shared memory corresponding 
to the edge adjacent to its parent p' in the spanning forest 
(any particle q can then easily check if p is a child of q ). 

Initially all system particles, except the seed, are inactive. 
In a nutshell, the particles that are touching the seed or 
other retired particle become roots; the root particles move 
around the partially constructed shape structure in a clock¬ 
wise manner until they find a valid position on the snake 
and become retired; follower particles follow the movement 
of the respective root until they become roots themselves. 
As we will see later, the initial snapshots of Figures [2] and [3] 
illustrate the spanning forest formation for the respective 
initial particle configurations. 

2.2 Hexagonal Shape Formation 

















Algorithm 1 Spanning Forest Algorithm for Shape Formation 


Depending on p’s current state, a particle p behaves as described below: 
inactive: If p is connected to a retired particle, then p becomes a root particle. Otherwise, if an adjacent particle p is a 
root or a follower, p sets the flag p.parent on the shared memory corresponding to the edge to p and becomes 
a follower. If none of the above applies, it remains inactive. 

follower: If p is contracted and connected a retired particle, then p becomes a root particle. Otherwise, it considers 
the following three cases: (i) if p is contracted and p’s parent p' is expanded, then p expands in the direction 
given by p.parent in a handover with p', and may need to adjust p.parent to still point to particle p after the 
handover; ( ii ) if p is expanded and has a contracted child particle p x , then p executes a handover with p ; (Hi) 
if p is expanded, has no children, and p has no inactive neighbor, then p contracts, 
root: Particle p runs the corresponding snake formation algorithm (Algorithm I or [3] for HEX or TRI resp.), and 

becomes retired accordingly. Otherwise, it considers the following three cases: (z) if p is contracted, it tries to 
expand in the direction given by RootDirection (p); (ii) If p is expanded and has a child p', then p executes 
a handover contraction with p ; (Hi) if p is expanded and has no children, and no inactive neighbor, then p 
contracts. 

retired: p performs no further action. 


RootDirection (p): 

Let i be the label of an edge connected to a retired particle, 
while edge i points to a retired particle do 

i label of next edge in clockwise direction 

return i 


We now investigate the Hexagonal Shape Formation (HEX) 
problem where the system of particles has to assume the 
shape of a hexagon (but for the outer layer, which may not 
be completely full) in G eq t ■ The hexagon will be constructed 
around the seed particle. Note that a hexagon in Geqt is ac¬ 
tually a disk, since it can be defined by the set of all nodes 
of G eq t within a certain distance r from a seed node. 

We will organize the particles according to a spiral snake 
structure which will incrementally add new layers to the 
hexagon, scaling naturally with the number of particles in 
the system. In order to characterize the snake formation for 
a given shape formation problem, one only needs to specify 
the direction in which the line of particles forming the snake 
should continue to grow, for each new particle added to the 
snake. Hence once a particle finds the next valid position on 
the snake, it will become retired and set the snake direction 
accordingly (by correctly setting the flag p.snakedir on the 
respective edge). Different rules for snake formation will 
realize different shapes. In particular, Algorithm [2] specifies 
the rules for the spiral snake formation for HEX. 

Initially, the seed particle p sets the flag p.snakedir in the 
shared memory corresponding to one of its adjacent edges 
(e.g., the edge with label 0). Any particle adjacent to a re¬ 
tired particle becomes a root following the spanning forest 
algorithm. Each root p moves in a clockwise fashion around 
the structure of retired particles until it finds the next po¬ 
sition to extend the hexagonal snake (i.e., a position con¬ 
nected to a retired particle via an edge flagged p.snakedir) 
and becomes retired, following Algorithm [5] (see Part (d) of 
Figure [l} . 

Figure [2] depicts some snapshots of a run of HEX algo¬ 
rithm. See Appendix for the proof of the following theorem. 

Theorem 1. Our algorithm solves the HEX problem in 
worst-case optimal 0(n 2 ) work. 

2.3 Triangular Shape Formation 

Now we investigate the Triangular Shape Formation prob¬ 
lem (TRI) where the system of particles has to assume a final 
triangular shape on G eq t (but for possibly the outer layer). 


Algorithm 2 Snake Formation for HEX 
if p is a contracted root then 

if p has an adjacent edge i to p with a flag p .snakedir, 
where p is retired then > retired condition 

while edge i is connected to a retired particle do 
i <— label of next edge in clockwise direction 
p sets the flag p.snakedir for edge i 
p becomes retired. 


As we discussed for the HEX problem, in order to solve 
the TRI problem in our Spanning Forest + Snake Forma¬ 
tion algorithmic framework, one only needs to setup the cor¬ 
rect rules for growing a ’’triangular snake”, which will be 
accomplished by Algorithm [3] The snake formation rules 
for the TRI problem are complex than the ones we had 
for the HEX problem, since we will need to explicitly take 
into account the formation of different layers of particles as 
we build the triangular structure (this was implicitly taken 
care by the spiral formation in the HEX algorithm). The 
TRI snake construction will start from the seed particle p, 
which will occupy one of the triangle corners. The seed will 
mark two of its adjacent edges as the direction along which 
two of the borders of the triangle will be formed, by setting 
p.border [left] and p.border [right] flags on the corresponding 
edges (we arbitrarily pick the edges will labels 0 and 1 out of 
p in our algorithm). These directions will be propagated by 
the particles that end up on one of the two sides. The seed 
starts the snake formation by setting the flag p.snakedir on 
its 0-labeled edge. From there on, Algorithm [3] will build 
the triangle snake layer by layer, alternating going ”to the 
left” and ”to the right”. Every time the snake touches one 
of the borders (Case 2 of Algorithm |3| , it sets up the rules 
for starting a new layer by setting the snake direction flags 
accordingly, first on the last particle of a layer (the one that 
just touched the border, in Case 2) and then on the first 
particle of the newly formed layer (Case 3). If a new layer 
is not needed, the snake proceeds to fill additional positions 
on the current layer (Case 1). Figure [3] illustrates this ap- 








proach through some snapshots of the execution of the TRI 
algorithm. The proof of the following theorem appears in 
the Appendix. 

Algorithm 3 Snake Formation for TRI 
if p is a contracted root then 

if p has an adjacent edge i to p with a flag p .snakedir, 
where p is retired then > retired condition 
bordertype — Border (p) 
if bordertype = null then 

> Case 1: continue on the same layer 
p sets p.snakedir for edge opposite to i 
(i.e., edge (i + 3) mod 6) 
else 

Let q be the border particle connected to p 
Let j be the edge of p opposite to the edge 
connecting p to q 

p sets p.border [bordertype] on edge j 

if p 7^ q then 

> Case 2: start a new layer 

p sets p.snakedir for edge j 

else 

> Case 3: snake direction from border 

if bordertype — left then 

p sets p.snakedir for edge (i + 5) mod 6 

else 

p sets p.snakedir for edge (i +1) mod 6 
p becomes retired 

Border (p): 

if p has an adjacent edge Ho a particle q with a flag 
q.border [bordertype], where bordertype £ {left,right} 

then 

return bordertype 

else 

return null 

Theorem 2. Our algorithm correctly solves the TRI prob¬ 
lem in worst-case optimal 0(n 2 ) work. 

3. CONCLUSION 

We presented a general algorithmic framework for shape 
formation problems in SOPS that combines our spanning 
forest algorithmic primitive with a snake formation prim¬ 
itive. We have shown that by carefully determining how 
to grow the appropriate snake structure, we were able to 
solve the HEX and TRI problems. We can easily to extend 
our snake primitive to build other shapes, such as square or 
rectangular shapes, or diamonds. It would be interesting to 
characterize all the general shapes that could be solved with 
our approach on G e q t (and possibly also for other infinite 
regular grid graphs, namely the square grid graph and the 
hexagonal grid graph, if we considered those as the under¬ 
lying graph G in the geometric amoebot model). Finally, 
we would like to evaluate the performance of our algorithms 
in terms of the worst-case number of asynchronous rounds 
necessary for termination. 
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Figure 2: Snapshots of the HEX algorithm. The seed is green, retired particles are black, roots are red and followers 
are blue. For a full simulation run of the algorithm see http://sops.cs.upb.de. 







Figure 3: Snapshots of the TRI algorithm. The seed is green, retired particles are black, roots are red and followers 
are blue. For a full simulation run of the algorithm see http://sops.cs.upb.de. 














APPENDIX 
A. ANALYSIS 

A. 1 Analysis of the Spanning Forest Algorithm 

We say followers and roots particles are active. As speci¬ 
fied in Algorithm^ only followers can set the flag p.parent. 

The first three lemmas demonstrate some properties that 
hold during the execution of the spanning forest procedure 
and will be used later to analyze the proposed algorithms 
for HEX and TRI problems. 

Lemma 1. For a follower p, the node indicated by p.parent 
is occupied by an active particle. 

Proof. Consider a follower p in any configuration during 
the execution of Algorithm [l] Note that p can only become 
a follower from an inactive state, and once it leaves the fol¬ 
lower state it will not switch to that state again. Consider 
the first configuration ci in which p is a follower. In the 
configuration Co immediately before ci, p must be inactive 
and it becomes a follower because of an active particle p 
occupying the position indicated by p.parent in co. The 
particle p is still adjacent to the edge flagged by p.parent in 
ci. Now assume that p.parent points to an active particle 
p' in a configuration c*, and that p is still a follower in the 
next configuration c^+i that results from executing an ac¬ 
tion a. If a affects p and p , the action must be a handover 
in which p updates its flag p.parent such that p.parent may 
be moved to the edge that now connects p to p in c*+i. If 
a affects p but not p 7 , it must be a contraction in which 
p.parent does not change and still points to p'. If a affects 
p but not p, there are multiple possibilities. The particle 
p might switch from follower to root state, or from root to 
retired state, or it might expand, none of which violates the 
lemma. Furthermore, p might contract. If p.parent points 
to the head of p 7 , p is still adjacent to the edge flagged by 
p.parent in Ci+i. Otherwise, p is a child adjacent to the tail 
of p in a and therefore the contraction must be part of a 
handover. As p is not involved in the action, the handover 
must be between p and a third active particle p". It is easy 
to see that after such a handover p.parent points to either 
p' or p". Finally, if a affects neither p nor p 7 , p.parent will 
still point to p 7 in q+i. □ 

Based on Lemma [l] we define a directed graph A(c) for 
a configuration c as follows. A(c) contains the same nodes 
as the nodes occupied in G eq t by the set of particles in c. 
For every expanded particle p in c, A(c) contains a directed 
edge from the tail to the head of p, and for every follower 
p in c, A(c) contains a directed edge from the head of p to 
p .parent. 

Lemma 2. The graph A(c) is a forest, and if there is at 
least one active particle, every connected component of in¬ 
active particles contains a particle that is connected to an 
active particle. 

Proof. In an initial configuration co, all particles are in¬ 
active and therefore the lemma holds trivially. Now assume 
that the lemma holds for a configuration a. We will show 
that it also holds for the next configuration c*+i that results 
from executing an action a. If a affects an inactive particle 
p, this particle either becomes a follower or a root. In the 
former case p joins an existing tree, and in the latter case 


p forms a new tree in A(q+i). In either case, A(c*+ 1 ) is 
a forest and the connected component of inactive particles 
that p belongs to in ci is either non-existent or connected to 
p in ci+ 1 . If a affects only a single particle p that is in state 
follower, this particle can contract or become a root. In the 
former case, p has no child p such that p .parent is the tail 
of p and also p has no inactive neighbors. Therefore, the 
contraction of p does not disconnect any follower or inactive 
particle and, accordingly, does not violate the lemma. In 
the latter case, p becomes a root of a tree which also does 
not violates the lemma. If a involves only a single particle 
p that is in state root, p can expand or contract or become 
retired. An expansion and becoming retired trivially cannot 
violate the lemma and the argument for the contraction is 
the same as for the contraction of a follower above. Finally, 
if a involves two active particles in a, these particles perform 
a handover. While such a handover can change the parent 
relation among the nodes, it cannot violate the lemma. □ 

The following lemma shows that the spanning forest al¬ 
ways makes progress, by showing that as long as the roots 
keep moving, the remaining particles will eventually follow. 

Lemma 3. An expanded particle eventually contracts. 

Proof. Consider an expanded particle p in a configura¬ 
tion c. Note that p must be active. If there is an enabled 
action that includes the contraction of p, that action will 
remain enabled until p eventually contracts when p is vali¬ 
dated in the current round. So assume that there is no en¬ 
abled action that includes the contraction of p. According to 
Lemma [2] and the transition rule from inactive to active par¬ 
ticles, at some point in time all particles in the system will 
be active. If the contraction of p becomes part of an enabled 
action before this happens, p will eventually contract. So as¬ 
sume that all particles are active but still p cannot contract. 
If p has no children, the isolated contraction of p is an en¬ 
abled action which contradicts our assumption. Therefore, 
p must have children. 

Furthermore, p must read at least one child p having its 
p' .parent flag pointing towards p over its tail and all chil¬ 
dren having their parent flags pointing towards p’s tail must 
be expanded as otherwise p could again contract as part of 
a handover. If p would contract, a handover between p and 
p would become an enabled action. We can apply the com¬ 
plete argument presented in this proof so far to p 7 and so on 
backwards along a branch in a tree in A(c) until we reach 
a particle that can contract. We will reach such a particle 
by Lemma |2] Therefore, we found a sequence of expanded 
particles that starts with p and ends with a particle that 
eventually contracts. The contraction of that last particle 
will allow the particle before it in the sequence to contract 
and so on. Finally, the contraction of p will become part 
of an enabled action and therefore p will eventually con¬ 
tract. □ 

A.2 Hexagonal Shape Formation Analysis 

Here, we show that the algorithmic primitives proposed 
in Section [2.2| solve the HEX problem correctly. 

Theorem 3. Our algorithm solves the HEX problem. 

Proof. We need to show that the algorithm terminates 
and that when it does, the system is in the shape of a 
hexagon. According to Lemma [2] every particle p eventu¬ 
ally activates. According to the spanning forest algorithm, 


if p is adjacent to the retired structure (initially the struc¬ 
ture only contains the seed particle), it becomes a root and 
moves in a clockwise manner around the retired structure 
until it eventually reaches the valid position that can extend 
the hexagon and becomes retired. By contradiction, assume 
p never becomes retired. Since the number of particles is 
bounded (and therefore the size of the formed retired struc¬ 
ture is bounded), there must be an infinite number of config¬ 
urations Ci where p had a root particle blocking its desired 
clockwise movement around the hexagona retired structure. 
Let p be the last root p sees as its clockwise neighbor over 
the retired structure (since once a particle becomes a root, 
it will stay connected to the hexagonal retired structure and 
always attempt to move in a clockwise manner, p is well- 
defined). Applying the same argument inductively to p ', we 
will get an infinite sequence of roots on the retired structure 
that never touch a valid spot pointed by q.snakedir flag of 
an already retired particle q , a contradiction, since the cur¬ 
rent retired structure (and the number of retired particles) 
is bounded. Therefore, every root eventually changes into a 
retired state. From Algorithm^ every follower in the neigh¬ 
borhood of a retired particle becomes a root. For every root 
q with at least one follower child, let c be the first configu¬ 
ration when q becomes retired. If q still has any child in c 
then all of its children p become roots. Applying this argu¬ 
ment recursively we will reach to a configuration such that 
there exists no root q having a follower child which proves 
that eventually every follower becomes a root. Putting it 
all together, eventually all particles become retired and the 
algorithm terminates. 

Note that it also follows from the argument above that 
the set of retired particles at the end of the algorithm forms 
a connected structure (since the particles start from a con¬ 
nected configuration and never get disconnected through the 
process). 

Now, we need to prove the correctness, i.e., that the re¬ 
sulting structure of retired particles is in a hexagonal form. 
Initially the hexagon only contains the seed particle, there¬ 
fore the claim holds trivially. By induction, let’s assume c 
is the first configuration in which the current formed struc¬ 
ture of the retired particles contains k retired particles and 
by induction hypothesis, assume that those particles form a 
valid hexagonal shape using k particles. According to Al¬ 
gorithm]^ the only way a root p can become the (k + l) th 
retired particle during or after c, is if it occupies the next 
valid position pointed by the flag q.snakedir , where q was 
the k -th particle to join the hexagonal shape. According 
to induction hypothesis, the k first retired particles form 
a hexagonal shape. By pointing to the next adjacent po¬ 
sition in counter-clockwise direction around the outermost 
retired particles in the current hexagonal structure, the flag 
q.snakedir points to the next position (according to counter¬ 
clockwise direction) on the last formed layer of the retired 
structure, or to a starting position on the next layer once 
the current layer is full, proving the correctness of the con¬ 
structed shape. □ 

We would like to measure the amount of the work of the 
proposed algorithm. 

Lemma 4. The worst-case work required by any algorithm 
to solve the HEX problem is 0(n 2 ). 

Proof. Consider a line of n particles on G eq t , where the 
seed particle is located on one end of the line, as an initial 


configuration of the particles. We label the particles con¬ 
nected to the seed starting with number 0 for the particle 
adjacent to the seed. The particle labeled i > 1 requires 
at least 2{i — 1 — \(i — 1)/Mj_i~|) > 2{i — 1 — |~(i — 1)/6~|) 
movements until it can lie contracted on the retired struc¬ 
ture where Mj 1 Mj > 6 and j > 1, indicates the capacity 
(i.e., the number of the retired particles) of the layer that the 
retired particle with label j belongs to. Therefore, any algo¬ 
rithm requires at least 2 Yl?=2 (* — 1 — |~(^ — 1)/6~|) = £4(n 2 ) 
work. □ 

Theorem 4. The algorithm proposed for HEX terminates 
in 0{n 2 ) work. 

Proof. To prove the upper bound, we simply show that 
every particle executes 0(n) movements. The theorem then 
follows. Consider a particle p. While p is in inactive or a 
retired state, it does not move. Let c be the first configura¬ 
tion when p becomes a follower. Consider the directed path 
in A(c) from the head of p to its root p . There always is 
such a path since every follower belongs to a tree in A(c) by 
Lemma |5] . Let P = (ao, ai,..., a m ) be that path in A(c) 
where ao is the head of p and a m is a child of p'. Accord¬ 
ing to Algorithm |T] p attempts to follow P by sequentially 
expanding into the nodes ao, ai,..., a m . The length of this 
path is bounded by 2 n and, therefore, the number of move¬ 
ments p executes while being a follower is 0(n). Once p 
becomes a root, it only performs expansions and contrac¬ 
tions around the retired structure until it reaches one of the 
valid positions on the hexagon. Since each root p and a 
retired particle q never connect from the same edge more 
than twice, and since the total number of retired particles is 
at most n, therefore the number of movements is bound to 
0(n) for p. Therefore, the number of movements a particle 
p totally executes is O(n), which concludes the theorem. □ 

A.3 Analysis for Triangular Shape Formation 

Now we need to show that the algorithmic primitives pre¬ 
sented in Section [T3| solve the TRI problem correctly. 

Theorem 5. Our algorithm solves the TRI problem. 

Proof. Again, we need to show that the algorithm ter¬ 
minates and that when it does, the system is the shape of 
a triangle. The termination part of the proof is identical 
to that for the HEX problem presented in Theorem [3] and 
hence it only remains to prove the correctness of the TRI al¬ 
gorithm. Assume we have three particles as the base case (to 
build the smallest size perfect triangle on Geqt)- The seed 
p* sets the p * .snakedir flag and the p * .border[left\ flag 
on its 0-labeled edge. A root particle q might have to move 
around the seed p* until it connects to edge 0 of the seed 
through an edges i. Since p sees both (border and snake) 
flags coming from the same particle, p becomes retired while 
it start constructing a new layer of the triangle and sets its 
p.snakedir flag such that the next particle continues filling 
this newly added layer (Case 3 of Algorithm [ 3 ]) . Particle p 
also sets p.border [left] appropriately to propagate the in¬ 
herited direction of the border from the seed to next layer. 
The only position that the third particle can stop on G eq t 
is the one pointed by p.snakedir and it is trivial to see that 
the resulting retired structure of the three particles is in a 
triangular shape. Let c be the first configuration in which 
the current formed structure of the retired particles contains 


k retired particles, and let q denote the ( k) th particle to be¬ 
come retired. By induction hypothesis, assume that those 
k particles form a triangle. According to Algorithm [ 3 ] the 
only way a root p can become the (k + l) th retired par¬ 
ticle during or after c, is if it occupies the valid position 
pointed by a flag q.snakedir. Depending on the location of 
q in the triangle, three cases may arise. First, consider the 
case when q is a left border particle (an analogous argument 
works if q is a right border particle). Since q is the last par¬ 
ticle added to the current valid triangular shape, we either 
have a perfect triangle after the addition of q or we have a 
perfect triangle plus particle q as the leftmost particle on 
a newly created layer. In the former case, given the next 
position pointed by q.snakedir , the root p follows Case 2 
of algorithm, which means that p will retire on the leftmost 
valid position on the next layer of the triangular structure, 
pointed by q.border[left]. In the latter, p follows Case 3 and 
will fill another position of the current layer next to q. In 
both cases the resulting retired structure still forms a valid 
triangular shape. Second, consider the situation where q is 
not a border particle (Case 1). Therefore, q is located on 
the last partially filled layer and q.snakedir is set to point to 
the next unoccupied snake spot on that layer, which is then 
filled by p, correctly extending the triangular structure, and 
proving the claim. □ 

Again, we would like to measure the work of the proposed 
algorithm. 

Lemma 5. The worst-case work required by any algorithm 
to solve the TRI problem is Q(n 2 ). 

Proof. With a very similar argument we had in Lemma|| 
one can verify that it is required to have at least 2 227=1 (i — 
1 — \(i — 1)/2~|) = f2(n 2 ) work for the algorithm to termi¬ 
nate. □ 

Theorem 6. The algorithm for TRI terminates in 0{n 2 ) 
work. 

Proof. Same argument we have in Lemma [4] holds here 
too. We just need to assume a triangular shape instead of a 
hexagonal one. □ 


